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We study disordered Peierls systems described by the Fluctuating Gap 
Model. We show that the typical electron states with energies lying deep 
inside the pseudogap are localized near large disorder fluctuations (instan- 
tons), which have the form of a soliton-antisoliton pair. Using the "saddle- 
point" method we obtain the average density of states and the average 

optical absorption coefficient at small energy. 
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Disorder in quasi-one-dimensional conductors strongly affects the electron states, re- 
ducing the tendency of these materials to develop 2k F - instabilities, such as the Peierls 
instability.0 Existing Peierls materials, like the conjugated polymer trans-polyacetylene, 
are known to suffer from various kinds of disorder: conformational defects, cross-links, 
impurities, etc.i While in a perfect Peierls chain the single-electron spectrum has a gap, 
with a value related to the amplitude of the periodic chain distortion, in a disordered 
chain the gap is filled. At weak disorder, the density of disorder-induced states that occur 
inside the gap is small, leading to a pseudogap, while at large disorder the gap disappears 
entirely and in the middle of the band the density may even diverge.! 

In this Letter, we find the typical form of the disorder-induced electron states lying 
close to the center of the pseudogap. In addition to giving new insight into these disor- 
dered systems, this result is useful, as it allows for a relatively easy calculation of disorder 
averages. As examples, we will explicitly calculate the average density of states and the 
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average absorption coefficient at small energies. 

It is crucial for our approach, that a large disorder fluctuation is required to create 
an electron state close to the center of the pseudogap, implying that the probability for 
such a state to occur is small at weak disorder. The main contribution to the density 
of states at small energy then comes from the disorder realizations close to one most 
probable fluctuation. The form of this fluctuation can be found consistently with the 
form of the wave function of the electron state induced inside the pseudogap. This wave 
function turns out to be localized in the vicinity of the disorder fluctuation. A similar 
situation is encountered when calculating the density of states of electrons moving in a 
random potential at large negative energiesBS These states were found to be localized in 
regions where the disorder potential has the form of a deep well. 

The approximate calculation of disorder averages, valid when the dominant contribu- 
tion comes from a small part of all possible disorder realizations, is closely related to the 
semiclassical approximation in quantum mechanics and field theory. In the path integral 
version of this approximation, the paths giving the largest contribution to the Green 
function lie close to the "saddle-point" of the Euclidean action, called the instanton.§ 

To describe the disordered Peierls chain we use the Fluctuating Gap Model (FGM).I 
It has previously been considered in the context of the thermodynamical properties of 
quasi-one-dimensional organic compounds (NMP-TCNQ, TTF-TCNQ),! and has been 
applied to study the effect of disorder on the Peierls transition,! as well as the effect 
of quantum lattice fluctuations on the optical spectrum of Peierls materials. In this 
model the electron motion is described by the one-dimensional Dirac equation, 



is the wave function of the single-electron state close to the Fermi energy ef = 0. The 





where o\ and 03 are the Pauli matrices and 




first term in the Hamiltonian h describes the free motion of the electrons and the two 
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amplitudes i/)r(x) and iPl(x) correspond to particles moving, respectively, to the right 
and to the left with the Fermi velocity Vp. The second term in the Hamiltonian describes 
the backward scattering of electrons from the lattice distortion wave, whose amplitude 
is proportional to A(x). Disorder is modelled by assuming A(x) to fluctuate randomly 
along the chain around some average value Ao, 

A(x) = A + V (x) , (2) 

where r](x) is the fluctuating part with a Gaussian correlator, 

(v(x)v(y)) = A5(x - y) . (3) 

In the absence of disorder (A(x) = A ), the electron spectrum has a gap between the 
energies e = —Aq and e = +Aq. At nonzero disorder the average density of states lying 
close to the middle of the pseudogap (\e\ <C A ) was found in Ref. |] by means of the 
"phase formalism" , 111 

(p(e)) <x \e\l-\ (4) 

where g = A /(vfA ) . We note that due to the charge conjugation symmetry (particle- 
hole symmetry) of the Dirac Hamiltonian h, the density of states is a symmetric function 
of the energy (p(—e)) = (p(e)), and in what follows we will assume e to be positive. 

From Eq.([|) it is clear that at weak disorder (j < 1) the density of states close 
to the middle of the pseudogap is strongly suppressed. As we mentioned above, the 
explanation for this is that a large fluctuation of A (a;) is required in order to create 
an electron state with energy £ < A . This motivates us to apply the "saddle-point" 
approach to study the typical electron states and to calculate disorder averages. The 
"saddle-point" disorder fluctuation (or instanton) fj(x) is the least suppressed one among 
the required large fluctuations. It can be found by minimizing, 

A [v(x)} = Jdxrf{x) - \i (e+ [v(x)\ - e) . (5) 

The first term in this equation describes the suppression of the probability of the fluc- 
tuation with the correlator Eq. (|]) (the weight p[i](x)] of the disorder configuration is 
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exp Jdxr] 2 (#))) , while the second term stems from the condition that the energy 
e + [t](x)} of the lowest positive energy single-electron state for the disorder realization 
r)(x) equals e. The factor fi is a Lagrange multiplier. The minimization of A [rj{x)} gives, 

fj(x) = fj,Ail)X{x)a 1 i() + (x) , (6) 

where ip+{x) is the wave function of the state with energy e + \f\{x)]. 

It may be shown by inspection that the solution of Eq. (§) is a soliton-antisoliton pair 
configuration, 



r)(x) = —VfK 



R \ { R 

tanh ( K(x — xq + — )J — tanh ( K{x — xq — — ) 



(7) 



where x describes the position of the disorder fluctuation in the chain, R is the distance 
between the soliton and the antisoliton, and K is determined by 

v F K = A tanh(fO?) . (8) 

The instanton is shown in Fig. 1 by plotting A(x) = Ao + fj(x). The spectrum of 
electron states that occur for this A(x) has previously been considered00 in relation with 
polarons in the Su-Schrieffer-Heeger model of conjugated polymers^S and is depicted 
in Fig. 2. It consists of a valence band (with highest energy — Ao), a conduction band 
(with lowest energy A ), and two localized intragap states with energies ±e + (R), where 

Thus, the soliton-antisoliton separation R is fixed by the condition £+(R) = £■ The 
two intragap states are the bonding and antibonding superpositions of the midgap states 
localized near the soliton and antisoliton: 



K 

4>±r(x) = ip*± L (x) = \ — 



± 



(10) 



cosh (k(x — xq — y)) cosh (K(x — Xq + ^)] 

(if)±(x)ip±(x) is schematically plotted in Fig.l). The energy splitting 2e decreases expo- 
nentially with the soliton-antisoliton separation, so that for e <^ Ao, 



J2«£ ln^, (11) 

where £ — Vp/A is the correlation length. The suppression factor Eq.(^), 

1 9 9 A 

A[v(x)]^—(2A fR^-\n^ , (12) 
2A g e 

also depends logarithmically on energy, so that the weight of the saddle-point configura- 
tion is, 
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p [t)(x)] oc 63 . (13) 

This result already gives a good estimate for the shape of the density of states inside the 
pseudogap at g -C 1 (cf. Eq.(|])). 

A more detailed calculation requires performing the Gaussian integration over the 
disorder realizations close to the "saddle-point" configuration Eq. ([?[). We then find the 
following expression for the average density of electron states per unit length, 

-l 
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For g<l, this agrees with Eq.(36) of Ref. [|, confirming the validity of the "saddle-point" 
approximation at small energies and weak disorder. The easiest way to get the result 



Eq.(|14|) is to use the correspondence between the averaging over disorder realizations r](x) 
and the quantum-mechanical averaging over the ground state for a certain double-well 
potential. The details of the calculation will be reported elsewhere. 

Having obtained the form of the most probable disorder-induced electron states, we 
can now also calculate in a relatively straightforward way the optical absorption coef- 
ficient for a half-filled chain at photon energy uj <C 2A and g <C 1. Again, only a 
large disorder fluctuation can make the energy difference between the empty and filled 
electron levels small. With the highest probability the photon absorption will induce a 
transition from the highest occupied to the lowest unoccupied electron state. Due to the 
particle-hole symmetry, the energy of the lowest unoccupied state at half-filling should 
equal +u/2 , while the energy of the highest doubly occupied state should equal —ui/2. 



Hence, the "saddle-point" disorder configuration, whose probability largely determines 
the absorption rate, is given by Eqs.(|7|), and (|) with e + (R) — u/2. 

Thus, in the "saddle-point" approximation, the absorption coefficient is the product 
of the averaged density of states (which is essentially the probability to find the necessary 
disorder fluctuation) and the strength of the optical transition between the two intragap 
levels: 

(a(u,)) = (C^E W» 6(E,-Eo-u>))m^- (p(e = w/2)> . (15) 

Here d is the electric dipole operator, |±) denote the wave functions of the intragap states 
(with energies ±a?/2) and C is an ^-independent coefficient (for small a>, we can neglect 
the weak cu-dependence of the real part of the dielectric constant). 

The wave functions of the intragap states |±) are the bonding and antibonding su- 
perpositions of the wave functions of the midgap states localized near the soliton and the 
antisoliton (see Eq.flTOD), from which the transition dipole matrix element is obtained as, 

(+\d\-) = q ~§ ■ (16) 

Here q denotes the electron charge and R is the soliton-antisoliton separation. Thus, for 
the asymptotic behavior of the averaged absorption coefficient at low photon energy we 
obtain, 

(a(u)) ocJ ( ln ^) 2 • ( 17 ) 

At this point we want to comment on the calculations of the optical conductivity 
in Refs. [5] and [11], in which the factorization approximation, (GG) = (G)(G), was used 



to evaluate the disorder average of the product of two Green functions. From above, 
it is clear that this approximation is not valid at low photon energies, as it results in 
the optical conductivity (as well as the absorption coefficient) being proportional to the 
second power of the weight Eq . (|13D , rather than the first (cf. Eq.(|T5|)). Of course, at 
weak disorder the absorption at photon energies uj 2A is small anyhow, but the 
factorization approximation makes it even much smaller. 



We would also like to point out that the small difference between the energies of the 
bonding and antibonding states is potentially dangerous for the saddle-point calculation 
of both the average optical absorption coefficient and density of states. The problem arises 
in the calculation of the contribution of the disorder realizations close to the "saddle- 
point" fluctuation: 

rj(x) = fj(x) + 5r)(x) . 

The perturbation Sh = aiSrj(x) can, in principle, strongly mix the bonding and antibond- 
ing states, because of the small energy denominator 2e appearing in the perturbation 
series. Such mixing would affect the values of both the energy splitting between the two 
states and the dipole matrix element. To see if that is the case, we considered the effective 
perturbation Hamiltonian, acting on the subspace of the two bonding and antibonding 
states, which includes the virtual excitations to all other (high-energy) electron states. 
We found that the off-diagonal matrix elements of this Hamiltonian are 0(e), which can- 
cels e in the denominator and makes the mixing of the two states small. This result is a 
direct consequence of the charge conjugation symmetry of the Dirac Hamiltonian Eq. ([]]). 
Thus, despite the small energy splitting, the saddle-point method is applicable. 

We conclude that in the FGM the most probable form of the wave function of the 
electron state lying deep within the pseudogap contains two peaks. The "saddle-point" 
disorder fluctuation, which induces such a state, has the form of a soliton-antisoliton pair 
and the peaks of the wave function are localized near the two kinks of this fluctuation (see 
Fig. 1). Away from the kinks, the electron wave function Eq.fllOl) falls off exponentially 
on a length scale £o- This observation is consistent with the fact, that the localization 
length at zero energy, 



hoc{e = 0) = 




calculated using the Thouless formula,li3 at weak disorder equals the correlation length 
£o- As we demonstrated, the instanton approach allows for a relatively easy calculation 
of the small-energy density of states and absorption coefficient. 
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Our results (Eqs.fll4"|) and Eq.(|I7|)) are valid if the density of the disorder-induced 
states is small, which is the case when |e| <C Aq and g ^ 1. It is useful, however, to 
comment briefly on effects of large disorder. For g ~ 1 the typical size of the disorder 
fluctuation i](x) on a scale of the correlation length £ becomes comparable to A , so 
that disorder fluctuations inducing the electron states with small energy are no longer 
suppressed, and for g > 2 the density of states even diverges at £ — 0.i This is essen- 
tially the singularity found long ago by DysonEl for a gapless system, because at strong 
disorder there is no principal difference between the electron states in Peierls insulators 
and conducting (A = 0) chains. In the latter case the localization length of the electron 
states diverges as e — > and the wave functions have an irregular structure, being large 
in many separated chain regionsJll Surprisingly, the "saddle-point" method gives the cor- 
rect exponent (2/g — 1) (cf. Eq. flTID ) for the energy dependence of the average density of 
states at all values of g. This suggests, that even at large g the typical form of the wave 



function may be close to the one given by Eq.(lC) in the regions where the wave function 
is large, and that multi-instanton disorder configurations (a gas of the soliton-antisoliton 
pairs) may become important at strong disorder. 
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FIGURE CAPTIONS 



FIG. 1. The form of A(x) = A + fj(x) for the instanton disorder fluctuation (thick 
line) and the electron density \tp + (x)\ 2 = \tp-(x)\ 2 for the corresponding intragap states 
(dotted line). 

FIG. 2. The spectrum of electron states for the instanton configuration A(x) plotted 
in Fig. 1. 
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